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We report the first measurement of the parity- violating asymmetry Apv in the elastic scattering 
of polarized electrons from 208 Pb. Apy is sensitive to the radius of the neutron distribution (R„). 
The result Apv — 0.656 ± 0.060 (stat) ± 0.014 (syst) ppm corresponds to a difference between 

the radii of the neutron and proton distributions R n — R p — 0.33 ^g'Jg fm and provides the first 

electroweak observation of the neutron skin which is expected in a heavy, neutron-rich nucleus. 

PACS numbers: 21.10.Gv, 21.65.Ef, 25.30.Bf, 27.80.-fw 190 < A < 219 



Nuclear charge densities have been accurately mea- 
sured with electron scattering and have become our pic- 
ture of the atomic nucleus, see for example In con- 
trast, our knowledge of neutron densities comes primar- 
ily from hadron scattering experiments involving, for ex- 
ample, pions protons or antiprotons @, 0], 
the interpretation of which requires a model-dependent 
description of the non-pcrturbative strong interaction. 
Due to the fact that the weak charge of the neutron is 
much larger than that of the proton, the measurement of 
parity violation in electron scattering provides a model- 
independent probe of neutron densities that is free from 
most strong- interaction uncertainties 0] . 

In the Born approximation, the parity violating cross- 
section asymmetry for longitudinally polarized electrons 
clastically scattered from an unpolarized nucleus, Apy, 
is proportional to the weak form factor Fy/(Q 2 ). This is 
the Fourier transform of the weak charge density, which is 
closely related to the neutron density, and therefore the 
neutron density can be extracted from an electro-weak 
measurement [al. 



A PV = 



°r-°l G F Q 2 F W (Q 2 



o-r + vl F ch {Q 2 ) 



(1) 



where <7r(l) is the differential cross section for elastic 
scattering of right (R) and left (L) handed longitudinally 
polarized electrons, Gf is the Fermi constant, a the fine 
structure constant, and F c h{Q 2 ) is the Fourier transform 
of the known charge density. However, the Born approx- 
imation is not valid for a heavy nucleus and Coulomb- 
distortion effects must be included. These have been ac- 
curately calculated [§] because the charge density is well 
known, and many other details relevant for a practical 
parity- violation experiment to measure neutron densities 
have been discussed in a previous publication (loj . 

One system of particular interest is the doubly-magic 
nucleus 208 Pb, which has 44 more neutrons than protons; 
some of these extra neutrons are expected to be found in 
the surface, where they form a neutron-rich skin. The 
thickness of this skin is sensitive to nuclear dynamics 
and provides fundamental nuclear structure information. 
A number of mean-field-theory models have been devel- 
oped that agree with the world's body of data on nuclear 



charge distributions and other nuclear properties lll-lla] . 
For 208 Pb, these are consistent with a radius of the point- 
neutron distribution R n between 0.0 - 0.4 fm larger than 
that of the point-proton distribution R p . In this paper 
we report a first measurement of Apy from 208 Pb, which 
is sensitive to the existence of the neutron skin. 

The value of the neutron radius of 208 Pb has important 
implications for models of nuclear structure and their ap- 
plication in atomic physics and astrophysics. There is a 
strong correlation between R n of 208 Pb and the pressure 
of neutron matter P at densities near 0.1 fm~ 3 (about 2/3 
of nuclear density) A larger P will push neutrons 

out against surface tension and increase R n . Therefore 
measuring R n constrains the equation of state (EOS), 
the pressure as a function of density, of neutron matter. 

The correlation between R n and the radius of a neutron 
star vns is also very interesting (l7| . In general, a larger 
R n implies a stiffer EOS, with a larger pressure, that 
will also suggest rjvs is larger. Recently there has been 
great progress in deducing rjvs from X-ray observations. 
From observations of X-ray bursts, Ozel et al. [3] find 
r/vs is very small, near 10 km, implying that the EOS 
softens at high density which is suggestive of a transition 
to an exotic phase of QCD. In contrast, Steiner et al. [lij 
conclude that r^s is near 12 km, leading to a prediction 
that Rn -Rp = 0A5± 0.02 fm for 208 Pb. This implies a 
stiffer EOS which leaves little room for softening due to 
a phase transition at high density. 

Recently Hebeler et al. [2(3| used chiral perturbation 
theory to calculate the EOS of neutron matter including 
important contributions from three-neutron forces. From 
their EOS, they predict R n - R p = 0.17 ± 0.03 fm for 
208 Pb. Monte Carlo calculations by Carlson et al. [2l[ 
also find sensitivity to three-neutron forces. The mea- 
surement of R n provides an important check of fun- 
damental neutron matter calculations, and constrains 
three- neutron forces. 

The EOS of neutron-rich matter is closely related to 
the symmetry energy S. There is a strong correlation be- 
tween R n and the density dependence of the symmetry 
energy dS/dp, with p as the baryon density. The symme- 
try energy can be probed in heavy-ion collisions |22j . For 
example, dS/dp has been extracted from isospin diffusion 
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data [23j using a transport model. 

The symmetry energy S helps determine the composi- 
tion of a neutron star. A large S at high density would 
imply a large proton fraction, which would allow the 
direct Urea process Q of rapid neutrino cooling. If 
Rn —Rp in 208 Pb were large, it is likely that massive neu- 
tron stars would cool quickly by direct Urea. In addition, 
the transition density from a solid neutron star crust to 
the liquid interior is strongly correlated with R n — R p 
25|. 

Rcinhard and Nazarewicz claim that R n — R p is tightly 
correlated with the dipole polarizability an [lil and 
Tamii et al. use this correlation to infer R n — R r , from a 
new measurement of an [2JJ . 

Atomic par ity violation (APV) is also sensitive to 

Rn [HHH 



29j . A future low-energy test of the standard 
model may involve the combination of a precise APV ex- 
periment along with PV electron scattering to constrain 
Rn [13 • Alternatively, measuring APV for a range of 
isotopes could provide information on neutron densities 

The measurement was carried out in Hall A at the 
Thomas Jefferson National Accelerator Facility. The ex- 
perimental configuration is similar to that used previ- 
ously for studies of the weak form factor of the proton 
and 4 He [3lTj33j . A 50 to 70 fiA continuous- wave beam 
of longitudinally polarized 1.06 GeV electrons was inci- 
dent on a 0.55 mm thick isotopically pure 208 Pb target 
foil. A 4 mm x 4 mm square beam raster prevented the 
target from melting. Two 150 /im diamond foils sand- 
wiched the lead foil to improve thermal conductance to 
a copper frame cooled to 20K with cryogenic helium. 
Elastically scattered electrons were focused onto thin 
quartz detectors in the twin High Resolution Spectrom- 
eters (HRS) [Hll. The addition of a pair of dipole sep- 
tum magnets between the target and the HRSs allowed 
us to achieve a forward scattering angle of 6\ a b ~ 5°. 
The HRS momentum resolution ensured that only elas- 
tic events (and a negligible fraction of inelastic events 
from the 2.6 MeV first excited state) were accepted by 
the quartz detectors. Cherenkov light from each quartz 
bar traversed air light guides and were detected by 2-inch 
quartz- window photo-multipliers (PMT). 

The polarized electron beam originated from a strained 
GaA sP p hotocathodc illuminated by circularly polarized 
light [35|. The accelerated beam was directed into Hall A, 
where its intensity, energy, polarization, and trajectory 
on target were inferred from the response of several mon- 
itoring devices. The sign of the laser circular polarization 
determined the electron helicity; this was held constant 
for periods of 8.33 ms, referred to as "windows" . The in- 
tegrated responses of detector PMTs and beam monitors 
were digitized by an 18-bit ADC and recorded for each 
window. Two "window quadruplet" patterns of helicity 
states (H h or — M — ) ensured complementary mea- 
surements at the same phase relative to the 60 Hz line 



power, thus canceling power-line noise from the asym- 
metry measurement. The right-left helicity asymmetry in 
the integrated detector response, normalized to the beam 
intensity, was computed for sets of complementary helic- 
ity windows in each quadruplet to form the raw asym- 
metry A raw . The sequence of these patterns was chosen 
with a pseudo-random number generator. 

Loose requirements were imposed on beam quality, re- 
moving periods of position, energy, or beam-intensity in- 
stability. No helicity-dependent cuts were applied, leav- 
ing a final data sample of 2 x 10 7 helicity- window quadru- 
plets. The design of the apparatus ensured that, after 
all corrections, the fluctuations in the fractional differ- 
ence of the PMT response between a pair of successive 
windows was dominated by scattered-electron counting 
statistics for rates up to 1 GHz. This facilitated the 
ability to achieve an Apy precision significantly better 
than 100 parts per billion (ppb) in a reasonable length 
of time. Careful attention to the desi gn a nd configura- 
tion of the photocathode laser optics [36| ensured that 
spurious beam-induced asymmetries were under control 
at this level. 

Random fluctuations in beam position and energy 
contributed the largest source of noise beyond count- 
ing statistics in A raw . Typical beam jitter in window- 
quadruplets was less than 2 parts per million (ppm) in 
energy, and 20 /im in position. This noise contribu- 
tion was reduced by measuring window differences Axi 
using beam position monitors and applying a correc- 
tion A^eam = ^CiAxi. The Cj's were measured several 
times each hour from calibration data in which the beam 
was modulated by using steering coils and an accelerat- 
ing cavity. The largest of the Cj's was ~ 50 ppm//im. 
The noise in the resulting A corr = A raw — Abeam was 
210 (180) ppm per quadruplet, for a beam current of 
50 (70) fj.A, dominated by counting statistics (~ 1 GHz at 
70 /LtA). Non- uniformities in target thickness due to ther- 
mal damage caused window-to- window luminosity fluctu- 
ations from variations in the target area sampled by the 
rastcred beam, leading to the degradation of A corr by 
~ 40%. This source of noise was eliminated by locking 
the raster pattern frequency to a multiple of the helicity 
frequency. Low-current calibration data, triggered on in- 
dividual scattered electrons, were regularly collected to 
evaluate the thickness of lead relative to diamond. 

Sensitivity of A corr to a transverse component of the 
beam polarization, coupled to the vector analyzing pow- 
ers (At) for 208 Pb and 12 C, was studied using special 
runs with fully transverse beam polarization. The sym- 
metry of the detector configuration as well as the mea- 
sured At values (to be published separately) resulted in 
an upper bound for a possible correction to A corr of 0.2%. 
The A raw and A corr window-pair distributions for the 
two complete data samples had negligible non-Gaussian 
tails over more than four orders of magnitude. To test 
the accuracy of error calculations and general statistical 
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behavior of the data, A corr averages and statistical er- 
rors were studied for typical one-hour runs, consisting of 
~ 50k quadruplets each. This set of 316 average A corr 
values, normalized by the corresponding statistical error, 
populated a Gaussian distribution of unit variance, as 
expected. 

A half- wave (A/2) plate was periodically inserted into 
the injector laser optical path, reversing the sign of the 
electron beam polarization relative to both the electronic 
hclicity control signals and the voltage applied to the po- 
larized source laser electro-optics. Roughly equal statis- 
tics were collected with this wavcplate inserted and re- 
tracted, suppressing many possible sources of systematic 
error. An independent method of hclicity reversal was 
feasible with a pair of Wien spin-rotators separated by a 
solenoid, providing an additional powerful check of sys- 
tematic control. Reversing the direction of the solenoidal 
field reversed the electron beam helicity while the beam 
optics, which depend on the square of the solenoidal 
magnetic field, were unchanged. The A/2 reversal was 
done about every 12 hours and the magnetic spin rever- 
sal was performed every few days. The dataset consisting 
of a period between two successive A/2 or magnetic spin- 
reversals is referred to as a "slug" . 

The spin reversals resulted in cumulative differences in 
beam position and energy of only 4 nm and 0.6 ppb re- 
spectively, leading to a run-averaged A aeam = — 39.0±5.9 
ppb. The asymmetry in beam charge, corrected by the 
intensity normalization of A raw , was 84. 0± 1.3 ppb, with 
the error determined using the correlation of measured 
beam intensity to PMT response which demonstrated 
the beam intensity monitors were linear to better than 
1.5%. Nonlincarity in the PMT response was limited 
to 1% in bench-tests that mimicked running conditions. 
As shown in Table HI the values of A corr are consistent 
within statistical errors for each of the reversal states. 
The reduced x 2 for A corr "slug" averages is close to 
one in every case, indicating that any residual beam- 
related systematic effects were small and randomized 
over the time period of A/2 reversals. The final result 
is A corr = 594 ± 50(stat) ± 9(syst) ppb where the sys- 
tematic uncertainty includes possible effects from A ueam , 
non-linearity in the detectors or beam charge monitors, 
and transverse asymmetry. The physics asymmetry Apy 
is formed from A corr by correcting for the beam polar- 
ization P u and background fractions fi with asymmetries 
At 



ipv 



i A corr — p u y) . Aifi 



Pb 



i-Ei/i 
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These corrections are summarized in Table [TT1 

The fraction of the accepted flux from 12 C in the de- 
tectors varied with time due to changes in the target; 
averaged over the run, the fraction / = (6.3 ± 0.6)%. 
The asymmetry of this background was determined to 
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TABLE I: Values of A corr and the statistical error, for each 
helicity reversal state and for the grand average. The \ 2 P er 
degree of freedom for each average is also shown. 



Correction 


Absolute (ppb) 


Relative(%) 


Beam Charge Normalization 
Beam Asymmetries A aeam 
Target Backing 
Detector Nonlinearity 
Transverse Asymmetry 
Polarization Pb 


-84.0 ± 1.5 
39.0 ± 7.2 
-8.8 ± 2.6 
± 7.6 
± 1.2 
70.9 ± 8.3 


-12.8 ± 0.2 
5.9 ± 1.1 
-1.3 ± 0.4 
± 1.2 
± 0.2 
10.8 ± 1.3 


Total 


17.1 ± 13.7 


2.6 ± 2.1% 



TABLE II: Corrections to Apy and systematic errors. 



be A C PV = 817 ±41 ppb using the Standard Model value 
for the e-N weak neutral isoscalar coupling and the mea- 
sured kinematics, with the uncertainty bounded by the 

4 He " 



precision measurement of Apy from 4 He [3l| . This was 
the only non-negligible background. An additional pos- 
sible systematic error in (Q 2 ) lay in the determination 
of the absolute value of 9i a b- A nuclear recoil technique 
using a water cell target (32| limited the scale error on 
(Q 2 ) to 1%. 

The spectrometer acceptance function e(8) character- 
izes the probability, as a function of scattering angle 9, 
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FIG. 1: Result of this experiment (red square) vs neutron 
point radius R„ in 208 Pb. Distorted-wave calculations for 
seven mean-field neutron densities are circles [38|], while the 
diamond marks the expectation for R n — R p . The blue 
squares show plane wave impulse approximation results. 
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for an electron to reach the detector after elastically scat- 
tering from 208 Pb. For example, the asymmetry averaged 
over the acceptance would be 



(A) 



Jd9 si 



sin 9 



j dO sin 



(3) 



where is the cross section. Sec Supplemental Ma- 
terial at http:/ /hallaweb.jlab.org/parity/prex/accept for 
the acceptance function e(6). The observed distribution 
of events corrected for the cross section, the background 
from the carbon (diamond) backing, and the effects of 
multiple scattering is used to extract e(9); corrections for 
energy loss in the target were negligible. To compare to 
predictions, one must integrate the theoretical asymme- 
try and the Q 2 over e(0). The systematic error in e(9) was 
evaluated from reasonable variations in the parameters of 
the simulation and resulted in an additional equivalent 
error in (Q 2 ) of 0.8%. Added in quadrature to the error 
arising from knowledge of (9) , we obtain an overall error 
in (Q 2 ) of 1.3%. We do not include this uncertainty in the 
total systematic uncertainty of the asymmetry. Using a 
calculation by Horowitz i, dA PV /dQ 2 is approximately 
30 ppm/GeV 2 , which would correspond to an additional 
systematic uncertainty on Apy of 3 ppb (0.5% of Apy). 

The beam polarization was continuously monitored by 
a Compton polarimctcr. Helicity-dependent asymme- 
tries in the integrated signal from backscattcred Comp- 
ton photons yielded P b = (88.2 ± 0.1 ± 1.0)% averaged 
over the duration of the run. The beam polarization was 
stable within systematic errors. An independent M0ller 
polarimctcr making nine measurements at different times 
during the run gave Pb = (90.3±0.1± 1.1)%. We used an 
average of these two measurements, Pb = (89.2 ± 1.0)% 
which conservatively accounts for the correlated system- 
atic errors between the two measurements. 

After all corrections, 
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Pb 
PV 



656 ± 60 (stat) ± 14 (syst) ppb 



at (Q 2 ) = 0.00880 ± 0.00011 GcV 2 . This result is dis- 
played in Figure [TJ in which models predicting the point- 
neutron radius illustrate the correlation of Apy and 
R n [HI- For each model, the calculation is performed 
using the neutron and proton weak charges q n = 0.9878 
and q p = —0.0721 and using the modeled neutron density 
but the experimental charge density. The importance of 
Coulomb distortions is emphasized by indicating results 
from plane- wave calculations, which are not all contained 
within the vertical axis range of the figure. A second- 
order polynomial fit over these models, as illustrated, 
implies a value for R n = 5.78 ^ ' 18 fm. Assuming a 
point-proton radius of 5.45 fm [37}, corresponding to the 
measured charge radius of 5.50 fm [l|, implies that the 
neutron distribution is 1.8er larger than that of the pro- 



future run is planned which will reduce the quoted un- 
certainty by a factor of three, to discriminate between 
models and allow predictions relevant for the description 
of neutron stars and parity violation in atomic systems. 
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0.33 
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fm [38| (see also [39(1). A 
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